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Abstract. A new family of parameters intended 
for composition studies is presented. They make 
exclusive use of surface data combining the in- 
formation from the total signal at each triggered 
detector and the array geometry. We perform an 
analytical study of these composition estimators in 
order to assess their reliability, stability and possible 
optimization. The influence of the different slopes of 
the proton and Iron lateral distribution function on 
the discrimination power of the estimators is also 
studied. Additionally, the stability of the parameter 
in face of a possible underestimation of the size of 
the muon component by the shower simulation codes, 
as it is suggested by experimental evidence, is also 
studied. 
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I. Introduction 

There are two main observation techniques of ultra- 
high energy cosmic rays, fluorescence and surface detec- 
tion, and they both have specific composition indicators. 
The most reliable technique at present for composition 
studies is fluorescence, where the longitudinal devel- 
opment of the charged component of the atmospheric 
shower is measured. Differences in composition, mani- 
fest themselves through differences in the cross section 
for interactions with atmospheric nuclei. These, in turn, 
are mapped as different depth of maximum development 
of the electromagnetic component in the atmosphere 
(X max ) and as dispersion in the position of this maxi- 
mum depth (AX max ) . If proton and Iron primaries are 
compared, the smaller cross section of the former will 
produce larger X max and AX max than for protons. 

On the other hand, surface detectors make a discrete 
sample of the shower front at ground level, measuring 
the shower particles density at different distances from 
shower core, which is called the lateral development 
of the shower. Beyond a few tens of meters from 
the shower axis, the particle content of the shower at 
ground level is dominated by just two components, 
electromagnetic (i.e., electrons, positrons and photons) 
and muonic. These two sets of particles propagate in 
a different way through the atmosphere: the electro- 
magnetic components propagates diffusively, while the 
muons do so radially from the last hadronic interaction 
that produce their parent mesons. Therefore, the main 
surface observables are related to the temporal and 



spatial distribution of particles. They are, for example, 
the slope of the lateral distribution function (LDF) used 
to fit the lateral development of the shower HI Q, the 
curvature of the shower front and several indicators of 
the time structure at a fixed point of the shower, like the 
rise time and fall time of the signal or their asymmetries 
0. Another promising parameter is the number of 
muons [4] (easily determined from scintillators but 
not from water Cherenkov tanks) because the differences 
between proton and Iron primaries are predicted to be 
around 70%. 

In general terms, fluorescence composition indicators 
are regarded as easier to observe and interpret, as well 
as less prone to systematic errors than surface param- 
eters do. However, fluorescence detectors suffer from a 
severely constrained duty cycle of approximately 10% 
of the total time available to surface detectors. This 
factor alone, which makes the statistics per unit time of 
surface arrays an order of magnitude larger than that of 
fluorescence detectors, gives a great attractive to search 
for reliable surface composition parameters, as we do in 
this work. 

Proton and Iron atmospheric showers are simulated 
using the AIRES Monte Carlo package (version 2.8.4a) 
[5| with QGSJET-II as the hadronic interaction model. 
We assume a surface array of of water Cherenkov 
detectors located on a triangular grid separated 1.5 km 
as in the Auger surface detector. The signal at each 
station is therefore the sum of the the electromagnetic 
and muonic components of the shower. The signal in 
each tank is measured in VEM (Vertical Equivalent 
Muon), i.e. the signal deposited by one vertical muon 
in a water Cherenkov tank, which is equivalent to 240 
MeV of deposited energy [6]. Therefore, to simulate the 
detector response easily we transform the energy of the 
shower particles at ground level into the measured signal 
considering that each muon contributes with 1 VEM and 
each electromagnetic particle of energy E, which are 
completely absorbed, with (E/240 MeV) VEM. 

II. Definition and optimization 

In the present work, we propose a new family of 
surface parameters defined as: 
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where the sum extends over all the triggered stations N, 
ro = 1000 m is a reference distance, Si is the signal in 
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VEM measured at the z — th station and Ti is the distance 
of this station to the shower axis in meters. 

First, we demonstrate that, under the present assump- 
tions, the primary identity discrimination power goes 
through a maximum around 6 = 3. Different detectors 
or a different layout could lead, in principle, to a 
different optimum value for b. The parameter Sb for a 
given event is constructed from the total signal in each 
triggered Cherenkov detector. Therefore, it depends on 
the normalization and shape of the lateral distribution 
function of the total signal. Close to the impact point 
of the shower, the signal is dominated by the elec- 
tromagnetic particles (photons, electrons and positrons) 
whereas at larger distances by muons. Fig. Q] shows the 
muon, electromagnetic and total signal in the Cherenkov 
detectors as a function of the distance to the shower 
axis for protons and iron nuclei. The zenith angle of the 
simulated events considered is such that 1 < sec (9 < 1.2 
and the primary energy 19 < log(_E/eV) < 19.1. Fig.Q] 
also shows the fits of the LDF of each component with 
a NKG-like function 01, 
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where we fix r s = 700 m and tq = 1000 m, and So and 
(3 are free fit parameters. If we consider proton and Iron 
primaries, the discrimination power of a mass sensitive 
parameter q, like Sb, can be estimated by using the so- 
called merit factor defined as, 
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where E[qA\ and l/ar^] are the mean value and 
the variance, respectively, of the distribution function 
of parameter qA with A = pr, fe. Assuming that the 
fluctuations of the total signal in an Auger Cherenkov 
detector are Gaussian, the distribution function for a 
given configuration of triggered stations is given by 
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where r; is the distance to the shower axis of the 
ith station (the first station, r\, is the closest one), 
S(r.i) is the average LDF evaluated at r"j, a[S(ri)] — 
1.06 [S^/VEM] 1 / 2 VEM |H| and f(r u . . . , r N ) is 
the distribution function of the distance of the differ- 
ent stations to the shower axis. Note that just two 
of the random variables {r±, ...,rjv} are indepen- 
dent, for instance, choosing to n and r2 (the first 
and second closest stations) as the independent ones, 
we can write f(n,...,r N ) = fia{ri,r 2 )S(r 3 - 
r3(n,r 2 )) ■ • ■ S(r N - r N (ri,r 2 )), where S(x) is the 
Dirac delta function. 

From Eqs. ([TJ and <j4j we obtain the expressions for 
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Fig. 1. Lateral distribution functions of the muon, electromagnetic 
and total signal in the Cherenkov detectors for simulated protons and 
Iron nuclei of 1 < sec0 < 1.2 and 19 < log(_B/eV) < 19.1. The 
hadronic interaction model used to generate the showers is QGSJET-II. 
The solid lines correspond to the fits with a NKG-like function (see 
Eq. (2j). 
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Here /i(rj) is the distribution function of the distance to 
the shower axis for the zth station and fi t j(ri, rj) is the 
distribution function of the distance to the shower axis 
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Fig. 2. i) as a function of b for vertical showers (1 < sec6 < 1.2) 
and 19 < \og(E/eV) < 19.1. Dashed and dotted lines correspond to 
the cases in which the slopes of the proton and iron LDFs are the same 
and equal to the corresponding to protons and iron nuclei, respectively. 




Fig. 3. Contour plot of rj(£, b)/r)(l, 3). 



of the ith and jth stations, 



conditions can be written as, 
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In order to simplify the expressions for the mean and 
variance of Sb we perform the following approximations, 



E[g( n )} - g(E[n}), 
cov[g(n),g(rj)] 
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where g{r) = S(r)(r/ro) b . 

We already have analytical expressions for the average 
LDFs of proton and Iron primaries obtained by fitting 
the simulated data, the other ingredients needed to 
calculate the mean value and the variance of Sb are the 
mean values of the distance to the shower axis for the 
different stations and the covariance between all pairs 
of those random variables. We obtain these quantities 
from a simple Monte Carlo simulation: we uniformly 
distribute impact points in a triangular grid of 1.5 km of 
spacing, like the Auger array, and then, for each event, of 
zenith angle such that sec 9 = 1.1 and azimuthal angle 
uniformly distributed in [0, 2ir], we calculate the distance 
of each station to the shower axis. Fig. |2] shows the 
discrimination power r\ as a function of b obtained under 
the mentioned assumptions and simplifications. We see 
that i) reaches the maximum at b = 3. 

III. Modifying the slope of the LDF 

We also study the discrimination power of Sb when 
the slope parameter (3 is modified but keeping constant 
the integrated signal for distances larger than the Moliere 
radius, tm — 80 m. The modified LDF that fulfills this 



S(r,(3) = 
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(r, + r M ) 
Beta(-r s /r M , -2(1 + 0), 1 + (3) Sp (r) 

where Beta(z, a, b) = f* dt t a (l-t) b and Sp (r) is the 
LDF of Eq. with the parameters Sq and /3o originally 
obtained from the fits. As a first calculation, fig. [2] also 
shows 7] as a function of b for the cases in which the 
slope parameter (3 of both LDFs is the same and equal 
to the corresponding to protons (dashed line) and Iron 
(dotted line). We see that r\ is considerable reduced when 
the slopes are equal. 

The slope of the proton LDF is smaller than the cor- 
responding to Iron (the absolute value is grater). Then, 
we modify the slope of both LDFs such that, (3 pr (£) = 
/3° r - (£- l)A/3 /2 and /3 /e (£) = p% + (£- 1) A/3 /2, 
where /3° r and /3® e are the proton and Iron slopes, 
respectively, obtained from the fits of the simulations, 
A/?o = P% - P a pr and £ is such that A/3(£) = £A/? , 
i.e., £ = 1 corresponds to the non modified case. Note 
that for £ = 0, fa - f3 fe = {(3 pr + p fe )/2. Fig. 
shows a contour plot of r](£,b)/r)(l,3) from where we 
see that as £ increases r] also increases. We also see 
that the maximum of r\ remains close to b = 3 almost 
independent of £. 

IV. Modifying the muon content of the 

SIMULATED SHOWERS 

There is experimental evidence of a deficit in the 
muon content of the simulated showers |9| . It is believed 
that such deficit is originated in the high energy hadronic 
interaction models which are extrapolations, over several 
orders of magnitude, of lower energy accelerator data. 
As mentioned, the total signal can be decomposed in the 
muon and electromagnetic signal. Therefore, to study 
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Fig. 4. Mean values of S3 for protons and iron nuclei as a function 
of / where / = 1 corresponds to the muon content predicted by 
QGSJET-II. 




how Sb changes as a function of the muon content of 
the showers we modify the total LDFs in the following 
way: S(r) — S em (r) + fS^(r) where / parametrizes the 
artificial variation in the muon component. Fig. [4] shows 
the mean values of S3 for protons and iron nuclei as 
a function of /. As expected, they increases with /. 
We also see that the iron curve increases faster than 
the proton one, which means that, for larger values of 
/, the discrimination power of S3 also increases. This 
happens because the muon content of the showers is very 
sensitive to the primary mass. Then, for large values 
of / the muon component becomes more important 
increasing the mass sensitivity of S3. 

Fig. [5] shows a contour plot of r)(f,b)/r](l,3) from 
which we see how the discrimination power of S& 
increases with the muon content of the showers and that 
the maximum is reached at b = 3 almost independently 
of /. 

V. Conclusions 

We propose a new family of parameters, which we call 
Sb, for composition analysis in cosmic ray surface de- 
tectors. The parameter is evaluated from the total signal 
and position of each triggered detector, on shower-to- 
shower basis. We perform an extensive analytical study 
of the most relevant properties of Sb- In particular, Sb has 
been optimized to distinguish between Iron and proton 
primaries assuming an Auger-like water Cherenkov de- 
tectors, showing that the discrimination power between 
both samples is maximum for 6 = 3. The potential dis- 
crimination power of the S3 parameter can be quantified 
by a merit factor 77 ~ 2.2 which is quite high if compared 
with the values attainable for competing parameters in 
current use. We have also demonstrated that, in case that 
the muon size is underestimated by simulation codes, as 
seems to be experimentally suggested, the parameter is 
not only stable but improves its discrimination power. 



Fig. 5. Contour plot of ri(f, b)/r)(l, 3). / = 1 corresponds to the 
muon content of the showers predicted by QGSJET-II. 
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